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Abstract. Wc suggest an alternative mathematical model for the electron in 
which the dynamical variables are a coframe (field of orthonormal bases) and a 
density. The electron mass and external electromagnetic field are incorporated 
into our model by means of a Kaluza-Klein extension. Our Lagrangian den- 
sity is proportional to axial torsion squared. The advantage of our approach is 
that it does not require the use of spinors, Pauli matrices or covariant differen- 
tiation. The only geometric concepts we use are those of a metric, differential 
form, wedge product and exterior derivative. We prove that in the special case 
with no dependence on the third spatial coordinate our model is equivalent to 
the Dirac equation. The crucial element of the proof is the observation that 
our Lagrangian admits a factorisation. 
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1. Introduction 

The Dirac equation is a system of four homogeneous linear complex partial dif- 
ferential equations for four complex unknowns. The unknowns (components of a 
bispinor) are functions of time and the three spatial coordinates. The Dirac equa- 
tion is the accepted mathematical model for an electron and its antiparticle, the 
positron, in a given external electromagnetic field. One of the main applications 
of the Dirac equation is spectral-theoretic: it determines with high accuracy the 
energy levels of the hydrogen atom. 

The geometric interpretation of the Dirac equation is rather complicated. It 
relies on the use of notions such as 

• spinor, 

• Pauli matrices, 

• covariant derivative (note that formula (|2.1j) for the covariant derivative of a 
spinor field is quite tricky). 
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There is also a logical problem with the Dirac equation in that distinguishing the 
electron from the positron forces one to resort to the concept of negative energy. 
Finally, the electromagnetic field is incorporated into the Dirac equation by means 
of a formal substitution which does not admit a simple geometric interpretation. 

The purpose of this paper is to formulate an alternative mathematical model 
for the electron and positron, a model which is geometrically much simpler. The 
advantage of our approach is that it does not require the use of spinors, Pauli 
matrices or covariant differentiation. The only geometric concepts we use are those 
of a 

• metric, 

• differential form, 

• wedge product, 

• exterior derivative. 

Our model overcomes the logical problem of distinguishing the electron from the 
positron: these correspond to clockwise and anticlockwise rotations of the coframe. 
And the electromagnetic field is incorporated into our model by means of a Kaluza- 
Klein extension which has a simple geometric interpretation. 

The paper has the following structure. In Section[2]we introduce our notation 
and in Section [3] we formulate the Dirac equation. In Section [4] we formulate our 
mathematical model and in Section [5] we translate our model into the language 
of bispinors. In Section [5] we prove Theorem 16.21 which is the main result of the 
paper: this theorem establishes that in the special case with no dependence on x 3 
our mathematical model is equivalent to the Dirac equation. The crucial element of 
the proof of Theorem 16. 21 is the observation that our Lagrangian admits a factori- 
sation; this factorisation is the subject of Lemma \6. II The concluding discussion 
is contained in Section [7l 



2. Notation and conventions 

Throughout this paper we work on a 4-manifold M equipped with prescribed 
Lorentzian metric g. All constructions presented in the paper are local so we do 
not make a priori assumptions on the geometric structure of spacetime {M, g}. 
The metric g is not necessarily the Minkowski metric. 

Our notation follows [HH]. In particular, in line with the traditions of particle 
physics, we use Greek letters to denote tensor (holonomic) indices. 

By V we denote the covariant derivative with respect to the Levi-Civita 
connection. It acts on a vector field and a spinor field as V a ii' 3 := d a v^ + r^ a7 w 7 
and 

V„£° := d a C + \a fi ad {d a a p bt + F 3 ^^)^ (2.1) 

respectively, where T 13 ai — {^j := \g l3S {d a g~fS + d 1 g a 5 — d$g ai ) are the Christof- 
fel symbols and ap are Pauli matrices. 



Dirac equation 



3 



We identify differential forms with covariant antisymmetric tensors. Given a 
pair of real covariant antisymmetric tensors P and Q of rank r we define their dot 
product as PQ := ^P ai ... ar Qi3 l ...p r g aiPl ■ ■■g arf3r . We also define ||P|| 2 := P P. 



3. The Dirac equation 

The following system of linear partial differential equations on M is known as the 
Dirac equation: 

<r aah (W + A) a% = me, a a ah (iV + A) a C = mVh- (3-1) 
Here £°, r]^ is a bispinor field which plays the role of dynamical variable (unknown 
quantity), m is the electron mass and A is the prescribed electromagnetic covector 
potential. The corresponding Lagrangian density is 



^Dir (£,??) 



Vaa aab V a% 



)-m(erj a +i\) V\detg\- (3.2) 



4. Our model 

A coframe i9 is a quartet of real covector fields j = 0,1,2,3, satisfying the 
constraint 

g = i9° (g> i9° - tf 1 <g) tf 1 - -& 1 ® -d 2 - i? 3 ® tf 3 . (4.1) 

For the sake of clarity we repeat formula (|4.1[) giving the tensor indices explicitly: 

n a — ?9° i9° — ?9 X ?9 X — ?9 2 ?9 2 — i9 3 ?9 3 

Formula (|4.ip means that the coframe is a field of orthonormal bases with 
orthonormality understood in the Lorentzian sense. Of course, at every point of 
the manifold M the choice of coframe is not unique: there are 6 real degrees of 
freedom in choosing the coframe and any pair of coframes is related by a Lorentz 
transformation. 

As dynamical variables in our model we choose a coframe d and a positive 
density p. These live in the original (1 + 3)-dimensional spacetime {M, g} and are 
functions of local coordinates (x 

In order to incorporate into our model mass and electromagnetic field we 
perform a Kaluza-Klein extension: we add an extra coordinate x and work on 
the resulting 5-manifold which we denote by M. We suppose that 

• the coordinate x A is fixed, 

• we allow only changes of coordinates (x°, x , x 2 , x 3 ) which do not depend 

4 

on x . 

We will use bold type for extended quantities. 
We extend our coframe as 

K = , *» = , (4-2) 
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(tf 1 + i$ 2 ) a = ^ + ^ 2 )^ e -2«- 4 7 (4.3) 

K = (° a ) (4.4) 

where the bold tensor index a runs through the values 0, 1, 2, 3, 4, whereas its 
non-bold counterpart a runs through the values 0, 1, 2, 3. In particular, the a in 
formula (|4.4[) stands for a column of four zeros. 

The coordinate x 4 parametrises a circle of radius ^ ■ Condition (|4.3[) means 
that the extended coframe # experiences a full turn in the (i? 1 , $ 2 )-plane as we 
move along this circle, coming back to the starting point. 

We extend our metric as 

/ 9af3 ~ ^2 A a Af) ^A a \ 

go* := • (4.5) 

V ^ -i / 

Formula (|4.5|) means that we view electromagnetism as a perturbation (shear) of 
the extended metric. Recall that in classical elasticity "shear" stands for "pertur- 
bation of the metric which does not change the volume" . It is easy to see that 
formula f|4. 5|) implies det g = — dct g , so dct g does not depend on A and, hence, 
the electromagnetic field does not change the volume form. 

Note that when i / the extended coframe and the extended metric no 
longer agree: 

g ^ ® i?° - i? 1 <8> tf 1 - ■d 2 ® tf 2 - tf 3 <g> tf 3 - i9 4 <g> i9 4 (4.6) 

(compare with (|4.ip ). The full physical implications of this discord are not dis- 
cussed in the current paper. We need the extended metric only for raising tensor 
indices (see formula (I4.9P below) and for this purpose the discord (|4. 6[) is irrelevant. 
We define the 3-form 

T ax := ~ (0° A dtf° - i? 1 A dti 1 — i9 2 A dti 2 - tf 3 A dtf 3 - fl 4 Arffl 4 J (4.7) 

=o 

where <i denotes the exterior derivative. This 3-form is called aa;«a/ torsion of 
the teleparallel connection. An explanation of the geometric meaning of the latter 
phrase as well as a detailed exposition of the application of torsion in field theory 
and the history of the subject can be found in [3]. For our purposes the 3-form (14.7[) 
is simply a measure of deformations generated by rotations of spacetime points. 
We choose our Lagrangian density to be 

m P ):=\\T™\\ 2 p (4.8) 

where 

||T ax ||2 ;= _L_ T a^ T a^ g « K g /3A gW (4 Q) 

Formula (14. 3|) implies 
i? 1 A dtf 1 + i9 2 A di9 2 = i? 1 A dtf 1 + i9 2 A cZi5» 2 - 4ml? 1 A i? 2 A i9 4 (4.10) 
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so our Lagrangian density p) does not depend on x 4 and can be viewed as 
a Lagrangian density in the original spacetime of dimension 1 + 3. This means, 
essentially, that we have performed a separation of variables in a nonlinear setting. 

Our action (variational functional) is J L(-d, p) dx° dx 1 dx 2 dx 3 . Our field equa- 
tions (Euler-Lagrange equations) are obtained by varying this action with respect 
to the coframe z9 and density p. Varying with respect to the density p is easy: this 
gives the field equation ||T ax || 2 = which is equivalent to L(-d,p) = 0. Varying 
with respect to the coframe $ is more difficult because we have to maintain the 
metric constraint (|4. 1[> ; recall that the metric is assumed to be prescribed (fixed). 

We do not write down the field equations for the Lagrangian density L{&, p) 
explicitly. We note only that they are highly nonlinear and do not appear to bear 
any resemblance to the linear Dirac equation (|3.ip . 



5. Choosing a common language 

In order to compare the two models described in Sections [3] and [4] we need to choose 
a common mathematical language. We choose the language of bispinors. Namely, 
we express the coframe and density via a bispinor field according to formulae 

S = eVa, (5.1) 



where 



p=\ s \^det\g a/3 \, (5.2) 

(ti + # 3 )« = \s\- 1 ev a J b , (5.3) 

(tf°-n* = N-y<w?^ (5-4) 

(^+i^) a = -\ s \-^a aab r l i > (5.5) 



Note that throughout this paper we assume that the density p does not vanish. 

Observe now that the right-hand sides of formulae (|5.2p - (|5.5[) are invariant 
under the change of bispinor field £ a i— » £ a e ltfi , rji i— > %e~ lv where <p : M — > R 
is an arbitrary scalar function. In other words, formulae (|5.2[) - (|5.5p do not feel 
the argument of the complex scalar s. Hence, when translating our model into the 
language of bispinors it is natural to impose the constraint 

lms = 0, s>0. (5.7) 



This constraint reflects the fact that our model has one real dynamical degree of 
freedom less than the Dirac model (seven real degrees of freedom instead of eight). 
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6. Special case with no dependence on x 3 

In addition to our usual assumptions (see beginning of Section [4]) we suppose that 

• the coordinate a; 3 is fixed, 

• we allow only changes of coordinates (x , x , x 2 ) which do not depend on a; 3 , 

• the metric does not depend on x 3 and has block structure 



/goo 


5oi 


502 


°\ 


5io 


5n 


512 





.920 


521 


522 













-v 



(6.1) 



the electromagnetic covector potential does not depend on a: 3 and has A3 = 0. 
We work with coframes such that 



tf 3 = 






W 



(6.2) 



'Zab 



1 1 2 



2 1 2 



(6.3) 

(6.4) 
This 

(6.5) 



We use Pauli matrices which do not depend on x 3 and take 

'1 

v -1 

We take 

% = £°*3oi ■ 

Then the scalar defined by formula (|5.1j) takes the form s 
scalar is automatically real and condition (|5.7p becomes 

l£T-l£T>o. 

It is easy to see that formulae (|6.1[) . (|6.3p - (|6.5p imply (|6.2j) . 

Formula (|6.4p means that our bispinor £ a , 77^ is determined by the spinor £ a . 
Thus, the spinor £ a becomes the (only) dynamical variable. We assume that this 
spinor does not depend on a: 3 . 

Observe that in the special case considered in this section both the Dirac 
model and our model have the same number of real dynamical degrees of freedom, 
namely, four. This is because under the assumption (|6.2p the coframe z9 and density 
p are equivalent to a spinor field £ a modulo sign (— £° gives the same & and p). 

Throughout this section summation is carried out either over indices 0, 1, 2 
or over indices 0, 1, 2, 4. In the latter case we use bold type. 

Put 



+ A a £ a a a ai e T mt a a 3ai) e\ V\^g\ ■ (6-6) 

The Lagrangian densities L^ iT (£) are formally related to the original Lagrangian 
density (|3.2p as follows: if we set rj^ = ±^°cr 3a j ) we get L mi (^, 77) = 2Lp ir (£). We say 
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"formally related" because in this section we assume that formula % = ±^ a cr 3Q j ) 
holds with upper sign, see 1|6.4|) . The Lj^ iT (£) and L^ ir (^) are, of course, the usual 
Dirac Lagrangian densities for an electron with spin up and spin down. 

Lemma 6.1. In the special case with no dependence on x 3 our Lagrangian den- 
sity factorises as 

Let us emphasise once again that throughout this paper we assume that 
the density p does not vanish. In the special case with no dependence on x 3 this 
assumption can be equivalently rewritten as 

^Dir(0 * WO (6-8) 

so the denominator in (|6.7[) is nonzero. 

Proof. Step 1. Let us show that it is sufficient to prove formula (|6.7|) under the 
assumption dA = 0, i.e. under the assumption that the electromagnetic covector 
potential A is pure gauge. Recall that dA stands for the exterior derivative of A. 

Suppose that we have already proved formula (|6.7p under the assumption 
dA = and are now looking at the case of general A. Let us fix an arbitrary 
point P on our 4-manifold M and prove formula (|6.7p at this point. To do this, 
we perturb the electromagnetic covector potential A in such a way that 

• A retains its value at the point P and 

• A satisfies the condition dA = in a neighbourhood of P. 

This can be achieved by, say, choosing some local coordinates on M and setting the 
components of A to be constant in this coordinate system. Now, this perturbation 
of the covector potential A does not change the LHS or the RHS of (|6.7|) at the 
point P because neither of them depends on derivatives of A. Hence, the case of 
general A has been reduced to the case dA = 0. 

Step 2. Let us show that it is sufficient to prove formula (|6.7p under the 
assumption A = 0. 

Suppose that we have already proved formula (|6.7j) under the assumption 
A = and are now looking at the case dA = 0. Let us modify the definition of the 
extended coframe by replacing (|4.4p with 

K = (~\ Aa ) ■ (6-9) 

In view of the condition dA — this modification of the extended coframe does not 
change axial torsion (|4.7|) but the extended coframe (|4.2p . (|4.3p . (|6.9[) now agrees 
with the extended metric (|4.5p : we have 

g = i?° ® i?° - tf 1 ® i? 1 - # 2 ® i? 2 - <8 i? 3 - i? 4 ® i? 4 (6.10) 
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as opposed to (J476J). Let us now perform a change of coordinates 

Z a =x a , a = 0,1,2,3, x i = x i -— A-dx. (6.11) 



?7i. 



Note that the integral / A ■ dx is (locally) well-defined because of the assumption 
dA = 0. The change of coordinates (|6.1ip is against the rules we stated in the be- 
ginning of Section Q] when describing our model (we changed the original Kaluza 
coordinate x 4 to a new coordinate x 4 ) but we are doing this only for the pur- 
pose of proving the lemma. In the new coordinate system i the extended coframc 
(|4.2p . (l4~3|) . ([63]) takes its original form ([4~2 |) -(j4~4 |) . the extended metric takes the 

form g a p = (compare with (|4.5[) ) and the electromagnetic covector 

V -v 

potential A is not affected (i.e. it has the same components in both coordinate 
systems). Observe now that in (|4.3[) we have retained the scalar factor e ~ 2lmx 
written in terms of the original Kaluza coordinate x A . Expressing x 4 in terms of 
x A in accordance with formula (|6. 11 [) we get 

(0 1 +l# 2 ) a = ^ e -2im^-2iSA-d x _ ( g 12) 

Let us now introduce a new coframc d in (l+3)-dimensional spacetime {M, g} 
related to the original coframe ■& as 

tf° = t?°, tf 3 = i? 3 , d 1 + id 2 = {d 1 + i$ 2 )e- 21 ! A ' dx . (6.13) 

Then formulae (|6"12l) . f6~Ll)l imply 

L0,p,O)=L(#,p;A). (6.14) 



Here L( ■ , • ; • ) is our Lagrangian density L( ■ , • ) defined by formulae 
(|4.7p ~ (|4.9p but with an extra entry after the semicolon for the electromagnetic 
covector potential. Formula (|6.14p means that in our model the introduction of an 
electromagnetic covector potential A satisfying the condition dA — is equivalent 
to a change of coframe (|6. 13|) . 

Formulae (|5.ip - (|5.6p . (|6.4p imply that the change of coframe (|6.13l) leads to 
a change of spinor field £ a = £ a e~ l f A-dx ^ Substituting the latter into (|6.6p we get 

^Dh.(e;0) = L± r (e;A). (6.15) 

Here Lp ir ( • ; • ) is the Dirac Lagrangian density L^ ir ( • ) defined by formula (|6.6p 
but with an extra entry after the semicolon for the electromagnetic covector po- 
tential. 

In the beginning of this part of the proof we assumed that we have already 
proved formula (|6.7p under the assumption A = so we have 

9 L Dir(£;°)- L Dir(£;°) 
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It remains to note that formulae (I6.14|) - (|6.16j) imply (|6.7j) . Hence, the case dA = 
has been reduced to the case A = 0. 

Step 3. In the remainder of the proof we assume that ^4 = 0. 

The proof of formula (|6.7p is performed by direct substitution: it is just a 
matter of expressing the coframe and density via the spinor using formulae (|5.1| - 
(|5.6p , (|6.4p and substituting these expressions into the LHS of (|6.7[) . However, 
even with A = this is a massive calculation. In order to overcome these technical 
difficulties we perform below a trick which makes the calculations much easier. 
This trick is a known one and was, for example, extensively used by A. Dimakis 
and F. Miiller-Hoissen 

Observe that when working with spinors we have the freedom in our choice of 
Pauli matrices: at every point of our (1 + 3)-dimensional spacetime {M, g} we can 
apply a proper Lorentz transformation to a given set of Pauli matrices to get a new 
set of Pauli matrices, with spinor fields transforming accordingly. It is sufficient to 
prove formula (|€>.7[) for one particular choice of Pauli matrices, hence it is natural 
to choose Pauli matrices in a way that makes calculations as simple as possible. 
We choose Pauli matrices 



aab 



jab 



# U S ■ 
"a *0a& 



'lab 



D 2ab 



i9 d s • 

u a *3a6 



(6.17) 



where 



'jab 



( S 0ab\ 
S lab 
S 2al> 

\ S 3ab/ 




(6.18) 



Here i9 is the coframe that appears in the LHS of formula (|6.7p . Let us stress that 
in the statement of the lemma Pauli matrices are not assumed to be related in any 
way to the coframe We are just choosing the particular Pauli matrices (|6.17l) . 
(|6.18p to simplify calculations in our proof. 

Examination of formulae (l5T|) - (|5^6l) . ([0]) . (163]) . (j6TT7jl . ([6T8]) shows that 
with our special choice of Pauli matrices we have £ 2 = whereas £ is nonzero 
and real. We are about to write down the Dirac Lagrangian density (|6.6p which is 



quadratic in £ so the sign of £ does not matter. So let £ a = 







where h : M 
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is a scalar function. We get 

j p 2h ■ 2h 

= —v a ai°p ad (d a ^ic + T\,a\ t ) + ... = —<j a ai ap a6 V a aP ld + ... 

j p 2h .... 

*H T^ct „ llvy „B _i_„a „ 12y7 „3 _j_„a „ 21v7 „8 ,„a _ 22vy _0 1 , 



[jaOa„ 11y7 „fl , o a0« 12*7 / 



12 



2/i 



+(<? 1 -t1? 2 ) 0, £r / 3 2l Va<^ ll +(i? 1 -** 3 ) a (r^ a4 Va<^ l4 ]+. 

8 

+ ( t 9 1 -^ 2 ) Q a0 2i V a i9 o/3 +(i9 1 - l ^ 2 ) Q ( T /3 22 V ct (i9 1 +^ 2 )^+. 
-(l? 1 -^^? 2 ) tt (^9 1 +^^? 2 ) /5 V Q ^+(^9 1 -^^9 2 ) a ^V ct (^9 1 +^^9 2 )' 3 ]+. 

+^ 2a ^V Q i9 0/3 +^ lQ ^V Q ^ 2 ' 3 -^ 2Q ^V Q i9 1/3 ]+. 

2h 

-^ 2Q 4V Q ^ 0/3 - t ? lQ ^V a ?9 2 ' 3 +79 2Q ^V Q ^ Vj ]+. 

= J [(t?° A d 1 ) ■ M 2 + (d 1 A d 2 ) • d-d + {d 2 A 0°) • dd 1 ] + ... 
8 

where the dots denote purely imaginary terms. Hence, 

2 (^o- tt oiVaS° " eV a a 6V a f ) = | A?? 1 ) • dtf 2 + (tf 1 A tf 2 ) • dtf° + (ti 2 A 0°) • dtf 1 ]. 
Formula (|6.6[) with A = can now be rewritten as 

- [(0° A tf 1 ) • dd 2 + (d 1 A i? 2 ) • + (tf 2 A 0°) • dz? 1 ] T m 



Put 



T ax := -(0° A - i? 1 A di? 1 - d 2 A d0 2 - 3 Arf0 3 ) 



p. (6.19) 
(6.20) 



(compare with (|4.7p ). The last term in (|6.20p vanishes in view of (|6.2[) . The coor- 
dinate x 3 is redundant so T ax can be viewed as a 3-form in (1 + 2)-dimcnsional 
Lorentzian space with local coordinates (x , x , x 2 ). Hence, we can define the scalar 

* T ax := i vTd^M (T ax )^£ Q/37 (6.21) 
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which is the Hodge dual of T ax . But y/\ detg\ e afjl = A 1 A i9 2 ) Q/ 3 7 so formula 
(|6.2ip can be rewritten as 

* T ax = T ax • (0° A 1 A tf 2 ) = -(0° A di}° - 1? 1 A di? 1 - d 2 A di? 2 ) • (0° A i? 1 A d 2 ) 

= \ [(0° A i? 1 ) • (M 2 + (tf 1 A tf 2 ) • d$° + (d 2 A • di} 1 ] . 
Substituting the latter into (|6.19[) we arrive at the compact formula 



_ * T ax T m 



(6.22) 



Substituting into the RHS of dHZZJ) we get 

32m LUfej 

9 "WO 

As our Lagrangian L($, p) is defined by formula 
been reduced to proving 



(* T ax ) 2 - — m 2 
V ' 9 



T ax || 2 = (*T ax ) 



ax\ 2 



16 



the proof of the lemma has 
(6.23) 



with ^4 = (recall that A initially appeared in the extended metric 
In view of P~2"]) . (f6T2")l formula (gTTJ) becomes 



(6.24) 



The difference between formulae (|6T2U)) and (|6~23|) is that the RHS of (|6~20"|) is 
expressed via the coframe § in the original spacetime whereas the RHS of (|6.24p 
is expressed via the coframe # in the extended spacetime, see (I4.2p ~ (l4.4p . In view 
of P~TU|) . (IQUf formula (f6T24)) can be rewritten as 

T ax = T ax + 1) At) A i? . (6.25) 

3 



The coordinate x 3 is redundant so T ax can be viewed as a 3-form in (1 + 3)- 
dimensional Lorentzian space with local coordinates (a; • X « X •} X 4 ). Hence, we 
can define the covector 

1 



(*T ax ) a :=- A AdiM(T ax ) a/37 £ a ^, 5 = 0,1,2,4, (6.26) 

which is the Hodge dual of T ax . It is easy to see that we have 

||T ax || 2 = -|| *T ax || 2 . (6.27) 

Note that in the LHS of (|6.27p we square a 3-form in (l+4)-dimensional Lorentzian 
space whereas in the RHS of (|6.27p we square a 1-form in (1 + 3)-dimensional 
Lorentzian space, so we took great care in getting the sign right. Substituting 
(|6~25|) into (jOB)) we get 

/4m o0\ 

(*T»),= (6.28) 
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where *T ax is the scalar defined by formula (|6.21|) . It remains to observe that 
formulae |Q7| . (f6T28|) imply ([Q3]) . □ 

The following theorem is the main result of our paper. 

Theorem 6.2. In the special case with no dependence on x 3 a coframe i? and a 
density p are a solution of the field equations for the Lagrangian density p) if 
and only if the corresponding spinor field is a solution of the field equation for the 
Lagrangian density L^ ir (^) or the field equation for the Lagrangian density L^^) . 

Proof. Denote by L(£) the Lagrangian density (|4.8|1 but with z? and p expressed 
via £. Accordingly, we rewrite the factorisation formula (|6.7ll as 

(0 ^"OF«)' (6 - 9) 

Observe also that the Dirac Lagrangian densities L^ ir defined by formula (|6.6p 
possess the property of scaling covariance: 

Ltir(e h = e 2h L± ir (t) (6-30) 
where h : M — > R is an arbitrary scalar function. 

We claim that the statement of the theorem follows from (|6.29[) and (|6.30l) . 
The proof presented below is an abstract one and does not depend on the physical 
nature of the dynamical variable £, the only requirement being that it is an element 
of a vector space so that scaling makes sense. 

Note that formulae (|6.29[) and (|6.30[) imply that the Lagrangian density L 
possesses the property of scaling covariance, so all three of our Lagrangian den- 
sities, L, Lp ir and LZ., have this property. Note also that if £ is a solution of 
the field equation for some Lagrangian density C possessing the property of scal- 
ing covariance then = 0. Indeed, let us perform a scaling variation of our 
dynamical variable 

£h-^ + /i£ (6.31) 

where h : M — ► R is an arbitrary "small" scalar function with compact support. 
Then = 5 J £(£) = 2 / h£(Q which holds for arbitrary h only if £(£) = 0. 

In the remainder of the proof the variations of £ are arbitrary and not nec- 
essarily of the scaling type (|6.3ip . 

Suppose that £ is a solution of the field equation for the Lagrangian density 
ijijjj.. [The case when £ is a solution of the field equation for the Lagrangian density 
ipjjj. is handled similarly] Then L^ ir (^) = and, in view of (|6.8|) . £fjir(£) ^ 0. 
Varying £, we get 

32m /" i 32m 



J^,(0 = 3 -^sJl^) 
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so 




32m 



9 




iir(e). 



(6.32) 



We assumed that £ is a solution of the field equation for the Lagrangian density 
L+ iT soSj L£ ir (£) = and formula fOgl implies that 8 J L(g) = 0. As the latter 
is true for an arbitrary variation of £ this means that £ is a solution of the field 
equation for the Lagrangian density L. 

Suppose that £ is a solution of the field equation for the Lagrangian density 
L. Then £(£) =0 and formula (j6T2"gi) implies that cither L+ ir (£) = or £ Dir (£) = 
0; note that in view of (|6.8p we cannot have simultaneously ip ir (£) = and 
-^Dir(£) = 0- Assume for definiteness that L^ iT (£) = 0. [The case when L^ iT (£) = 
is handled similarly] Varying £ and repeating the argument from the previous 
paragraph we arrive at (|6.32[) . We assumed that £ is a solution of the field equation 
for the Lagrangian density L so 6 f L(£) — and formula 1|6.32[) implies that 
6 J = 0. As the latter is true for an arbitrary variation of £ this means that 

£ is a solution of the field equation for the Lagrangian density L^ iv . □ 

The proof of Theorem 16.21 presented above may appear to be non-rigorous 
but it can be easily recast in terms of explicitly written field equations. 

7. Discussion 

The mathematical model formulated in Section 2] is based on the idea that every 
point of spacetime can rotate and that rotations of different points are totally inde- 
pendent. The idea of studying such continua belongs to the Cosserat brothers [7]. 
Recall that in classical elasticity the deformation of a continuum is described by a 
(co)vector function u, the field of displacements, which is the dynamical variable 
(unknown quantity) in the system of equations. Displacements, of course, gener- 
ate rotations: the infinitesimal rotation caused by a displacement field u is du, 
the exterior derivative of u. The Cosserat brothers' idea was to make rotations to- 
tally independent of displacements, so that the coframe (field of orthonormal bases 
attached to points of the continuum) becomes an additional dynamical variable. 

Our model is a special case of Cosserat elasticity in that we model spacetime 
as a continuum which cannot experience displacements, only rotations. The idea 
of studying such continua is also not new: it lies at the heart of the theory of 
teleparallelism (= absolute parallelism), a subject promoted in the end of the 
1920s by A. Einstein and E. Cartan [HHHl HO] - It is interesting that Einstein pursued 
this activity precisely with the aim of modelling the electron, but, unfortunately, 
without success. 

The differences between our mathematical model formulated in Section[4]and 
mathematical models commonly used in teleparallelism are as follows. 

• We assume the metric to be prescribed (fixed) whereas in teleparallelism it 
is traditional to view the metric as a dynamical variable. In other words, in 
works on teleparallelism it is customary to view (|4.ip not as a constraint but 
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as a definition of the metric and, consequently, to vary the coframe without 
any constraints at all. This is not surprising as most, if not all, authors who 
contributed to teleparallelism came to the subject from General Relativity. 

• We choose a very particular Lagrangian density (|4.8p containing only one 
irreducible piece of torsion (axial) whereas in teleparallelism it is traditional 
to choose a more general Lagrangian containing all three pieces (tensor, trace 
and axial), see formula (26) in [3]. 

We now explain the motivation behind our choice of the Lagrangian den- 
sity (|4.8[) . Suppose for simplicity that we don't have electromagnetism, i.e. that 
A = 0, in which case the extended coframe and extended metric agree (|6.10l) . Let us 
perform a conformal rescaling of the extended coframe: i?- 7 i— > e^i?- 5 , j = 0, 1, 2, 3, 4, 
where h : M — > R is an arbitrary scalar function. Then the metric and axial torsion 
scale as g i— » e 2h g and 

T ax | — ^ e 2h T ax ^ ^ 

respectively. Here the remarkable fact is that the derivatives of h do not appear 
in formula (|7.ip which means that axial torsion is the irreducible piece of torsion 
which is conformally covariant. It remains to note that if we scale the density p as 
p i— >■ e 2h p then the Lagrangian density (14. 8|) will not change. 

Thus, the guiding principle in our choice of the Lagrangian density (|4. 8|) is 
conformal invariance. This does not, however, mean that our mathematical model 
formulated in Section |4] is conformally invariant: formula (|4.4p does not allow for 
conformal rescalings. The Kaluza-Klein extension is a procedure which breaks 
conformal invariance, as one would expect when introducing mass. 

The main result of our paper is Theorem 16.21 which establishes that in the 
special case with no dependence on x 3 our mathematical model is equivalent to 
the Dirac equation. This special case is known in literature as the Dirac equation 
in dimension 1 + 2 and is in itself the subject of extensive research. 

This leaves us with the question what can be said about the general case, 
when there is dependence on all spacetime coordinates (x 3 ). In the gen- 

eral case our model is clearly not equivalent to the Dirac equation because it lacks 
one real dynamical degree of freedom, see last paragraph in Section[5] Our plan for 
the future is to examine how much our model differs from the Dirac model in the 
general case. We plan to compare the two models by calculating energy spectra of 
the electron in a given stationary electromagnetic field, starting with the case of 
the Coulomb potential (hydrogen atom). 

The spectral-theoretic analysis of our model will, however, pose a monumental 
analytic challenge. There are several fundamental issues that have to be dealt with. 

• Our model does not appear to fit into the standard scheme of strongly hy- 
perbolic systems of partial differential equations. 

• The eigenvalue (— bound state) problem for our model is nonlinear. 

• Our construction relies on the density p being strictly positive. This assump- 
tion may fail when one seeks bound states other than the ground state. 
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